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CONVECTIVE HEAT -TRANSFER COEFFICIENTS FROM A SOLUTION OF 
THE CONDUCTION EQUATION FOR A WALL SEPARATING TWO 
FLUIDS, ONE HAVING AN OSCILLATING TEMPERATURE 
by Ronald G. Huff 
Lewis Research Center 

SUMMARY 

The temperature response of a wall which separates two fluids, one having a sinus- 
oidally oscillating temperature has been analytically determined. Two mathematical 
solutions are presented; one uses thermal properties distributed through the wall; the 
other assumes that the wall thermal properties can be lumped. 

The solutions show that the convective heat -transfer coefficients on both sides of the 
wall can be determined experimentally if a single wall, the hot fluid, and coolant tempera- 
tures are known as a function of time. From traces of these temperatures the frequency 
of oscillation and either amplitude ratio or the phase lag angle between the forced fluid 
and wall temperatures are measured. These values are used in the solutions for the 
coefficients. The method is particularly applicable to rocket engines and is generally 
applicable to heat exchangers. 

Charts are included to allow graphical solutions for the distributed -wall -property 
case. Comparisons of the two mathematical solutions are made which indicate that it 
may be possible to use the simpler lumped-wall-property solution. This requires, among 
other thirds, thin walls and low forcing frequencies. The lumped -wall -property solution 
may be applied with greater accuracy when the wall temperature is measured at the mid- 
point of the wall. 


INTRODUCTION 

Transient and steady -state analyses have been used to design calorimeters for de- 
termining convective heat transfer in heat exchangers, rocket engines, and aerodynamic 
heat -transfer studies. The steady-state calorimeter of reference 1 makes use of the 
temperature gradient in a material of known conductivity and geometry; the transient 



type (ref, 2) uses the temperature response of a material to a change in driving tem- 
perature. 

The temperature response of a wall having one side insulated and the other exposed 
to a fluid with a sinusoidally varying temperature has been investigated by Anderson 
(ref. 3) and this author (ref. 4). Equations for the convective heat -transfer coefficients 
were derived in terms of the wall material properties, frequency, and phase lag angle 
between the. force fluid temperature and the responding wall temperature. The applica- 
tion of these equations is limited to walls having either one side insulated or having fluid 
temperatures and convective heat-transfer coefficients that are equal on both sides of the 
wall. 

Of more general interest is the problem of determining the heat-transfer coefficients 
in a heat exchanger or cooled rocket engine. This problem can also be solved, by sinus- 
oidally varying the temperature of the hot fluid or the coolant and measuring the tempera- 
ture response of the wall. The combustion temperature in a rocket engine, for example, 
can be oscillated by changing the ratio of the oxidizer -to -fuel mass flow rates while 
maintaining the total propellant mass -flow rate constant. This will require separate 
coolant and propellant systems so that the mass flow rates that affect the heat-transfer 
coefficients can be held constant. This uncoupled system is completely reasonable from 
a research point of view. 

The purpose of this analytical study was to derive the equations for the convective 
heat-transfer coefficients in terms of the sinusoidally varying temperature response of 
a wall which separates two fluids, one having a sinusoidally oscillating temperature. 

This study was performed at NASA Lewis in conjunction with a rocket nozzle heat -transfer 
program. 


STATEMENT OF THE PROBLEM 

Figure 1(a) shows a wall that separates two moving fluids. The hot-gas-side fluid 
temperature Tq is greater than the coolant temperature T^, causing heat to flow 
through the wall in the positive x direction. The problem is to determine the convective 
heat -transfer coefficients if either the hot -gas or coolant temperature is varied sinus- 
oidally and if the temperature response of the wall is measured at only one point. Fig- 
ure 1(b) shows what the wall temperature might look like at any given instant in time. 


METHOD OF ATTACK 

If the hot-gas temperature is made to vary sinusoidally, the wall temperature will 
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(b) Transient temperature response of plate to sinusoidally forced fluid temperature. 

Figure 1. - Continued. 

respond sinusoidally but will lag the driven temperature by an angle cp. In addition, the 
amplitude of the wall temperature oscillation will be less than that of the driven tempera- 
ture. Both the phase lag cp and the ratio of the amplitude of the wall temperature to 
the amplitude of the driven temperature 0/AT are, among other things, functions of the 
convective heat -transfer coefficients h. Finding the relation between the convective heat- 
transfer coefficients and the phase lag cp, or the amplitude ratio is required in 

order to solve for the coefficients. 

To find this relation, the temperature response of the wall to a sinusoidally driven 
fluid temperature is derived. From this solution and steady-state heat-transfer conditions 
(which relate the ratio of the convective heat-transfer coefficients to the ratio of the 
temperature drops between fluids and wall), the absolute values of the coefficients can be 
determined as functions of either phase lag or amplitude ratio. The quantities that 
must be known as a function of time are the hot-gas temperature Tq, the wall tempera- 
ture at any point x (T(x, t)) and the coolant temperature T^.. From these quantities 
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either the phase lag or the amplitude ratio can be determined as well as the ratio of the 
convective heat -transfer coefficients. With the ratio of the coefficients, the frequency of 
the temperature oscillation, the wall material properties, and either <P or S^/ATq, 
the convective heat -transfer coefficients h^ and h^ can be calculated. 

It is also possible to calculate the coefficients by oscillating the coolant temperature 
sinusoidally. This approach requires the measurement of the same quantities as when 
the hot fluid temperature is oscillated. 


FORMULATIONS AND SOLUTIONS OF THE PROBLEM 

Two analytical solutions have been found, each of which relate the heat -transfer 
coefficients to either the phase lag or the amplitude ratio. The first solution (case 1) uses 
the second-order partial differential equation for transient, one -dimensional, heat con- 
duction in a plate. The boundary conditions are that one surface is heated convectively 
and the other surface is cooled convectively. The second solution (case 2) uses a lumped 
wall property approach. The differential equation used in case 2 is obtained by equating 
the heat stored in the wall as a function of time to the difference between the heat trans- 
ferred to one side of the plate convectively and the heat taken from the other side of plate 
convectively. The case 2 solution for the wall temperature can be broken into three 
parts: the starting transient, an offset of the average value from the value before the 
start of the oscillation, and the steady-state oscillation. The steady-state oscillating part 
of the case 2 solution is used extensively in this work. 

Solutions for both case 1 and case 2 are given in this report so that the reader can 
determine for himself when the simpler, though less accurate, case 2 solution is applica- 
ble. 


Distribution Wall Properties (case 1) 

The following assumptions have been used to obtain the case 1 solution: 

(1) The heat flows through the plate in the x direction only (one -dimensional heat 
conduction) . 

(2) The wall or plate properties (density, specific heat, and thermal conductivity) 
are constant. 

(3) The convective heat-transfer coefficients on both sides of the wall are constant. 

(4) The coolant temperature is constant. 

Determining the validity of assumption (1) is left to the experimenter with his parti- 
cular application since no general comments are possible as to when one-dimensional 
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heat flow is present. The second assumption, that of constant wall properties, will de- 
pend on the wall material and the temperature through the wall including the amplitude of 
the oscillations. The wall temperature amplitude can be minimized so that the assumption 
of constant wall properties as a function of time should be applicable. If the wall is thin 
enough or the heat flux low enough, the wall properties may not change appreciably from 
one surface to the other. The third assumption, constant coefficients, appears justified 
because the coefficients are usually not strong functions of pressure or temperature. 

They are, however, strong functions of the mass velocity and for that reason the mass 
flow rate of the fluids are held constant. In the experimental rocket engine this is ac- 
complished by uncoupling the coolant flow rate from the propellant flow rate by use of 
separate systems. This allows a constant coolant flow rate to be set. The total mass 
flow rate of the propellants (oxidizer plus fuel) is held constant but in order to vary the 
combustion temperature the ratio of the individual propellant mass flow rate is varied. 

The validity of assumption (4) can be insured by limiting the amplitude of the hot fluid 
temperature oscillation to a value that forces the amplitude of the wall temperature os- 
cillation to be just adequately measureable. This measurement can be aided by the use 
of amplifiers that would allow the use of smaller hot -fluid amplitudes. With small wall- 
temperature amplitudes the coolant temperature will remain practically constant. 

The controlling differential equations and boundary conditions are given next. The 
one -dimensional transient heat conduction equation 

( 1 ) 

0x2 « 9'^ 

is solved for the temperature distribution in a wall which has convective heat transfer 
over its two surfaces (fig. 1). (Symbols are defined in appendix A. ) The hot-gas side 
fluid temperature is driven sinusoidally and is given in equation form by 

Tq = AT^e'^^^"^ + Tq (2) 

where Tq is given, in the rocket ei^ine application by the usual chemical calculational 
techniques and the time dependent term is sensed by the chamber pressure. It is not 
necessary to know the numerical value of ATq if the phase lag is used for calculating 
the coefficients since the phase determination is independent of the absolute value of the 
temperature oscillation. If the amplitude ratio is used for calculatit^ the coefficients 
the ATq must be computed from the change in chamber pressure. The boundary condi- 
tion at the hot -gas side of the wall surface x/L = 0 equates the heat transferred by con- 
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vection to that conducted away from the surface into the wall. This is given in equation 
form as 


to ■ (3) 

3x 

The boundary condition at the coolant surface (x/L = 1) equates the heat conducted to the 
surface to the heat transferred convectively to the coolant. In equation form this is given 
by 


h fT(L, t) - T 1 = - (4) 

The solution for the wall temperature as a function of time and location within the 

wall consists of the sum of the transient 0(x, t) and steady -state T (x) solutions. The 

s 

steady -state solution is found by settir^ 3T/3 t equal to zero in equation (1) and assum- 
ing the wall temperature at x/L = 0 and 1 to be given as and respectively. 

This yields the usual linear temperature distribution with distance and can be written as 

Tg(x) = (T^^ - Tq^) — + Tq^ (5) 

The transient solution e(x, t) is found by assuming the usual product solution: one 
a function of time only F(r), and the other a function of the location within the wall 
measured from the heated side of the wall X(x). Equations (1), (3), and (4) are modified 
to account for the change in variable from T(x, t) to 0(x, r). The details of the solution 
of equation (1) using its pertinent boundary conditions (eqs. (2) to (4)) are shown in 
appendix B. 

The transient part of the solution 6{x, t) is given in appendix B in complex form by 
equation (B25) presented here in functional form as 

e(x, t) = P (x, L, T], xjy^, )exp 

ATq 


where 


^ (l ’ 


1 



From this equation, the phase lag <p (eq. (B26)) and the amplitude of the wall tempera- 
ture oscillation (eq. (B27)) are derived. Equations (B28) to (B33) define the parameters 
used in equations (B25) to (B27). Equations (B25) to (B27) are rewritten in nondimensional 
form as equations (B34), (B36), and (B37), respectively. The governing nondimensional 
parameters which appear in the transient solution equation (B34) are xj^Q and the 
convective heat-transfer parameters; ?jL, the wall properties-thickness-frequency para- 
meter; and x/L, the location at which the wall temperature is being considered. The 
equations for phase lag (eq. (B36)) and amplitude ratio (eq. (B37)) are given as functions 
of the nondimensional parameters. Either one of these equations gives a relation between 
the desired coefficients but, since these equations are derived from equation (B34), it 
is impossible to solve them simultaneously for the desired coefficients. Instead another 
independent relation must be found. 

The required relation is derived from the steady-state (or mean) condition that re- 
quires the heat transferred convectively to the wall to equal the heat transferred con- 
vectively to the coolant. This relation is expressed mathematically as 


t^c _ '^G - 

T(L) - Tp 


( 6 ) 


From the definition of ;// 


h = (7) 

The ratio of convective heat -transfer coefficients in equation (6) is defined as R. Calcu- 
lating the ratio R using equation (7) gives 



(8) 


where R is also given by 
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( 9 ) 


- T(0) 

^ 

T(L) - 

(Eqs. (8) and (9) are also eqs. (B50) and (B49), respectively.) Since Tq, T^, and a 
wall temperature have been given which for the present is assumed to be constant allow- 
ing T(0) = T(L), R may be calculated from equation (9). With the value of R known, 
equation (8) then gives the second required relation between and Only one 

unknown exists in each of equations (B34), (B36), and (B37). Thus, and yp ^ can 
pbw be calculated either by knowing the phase lag (p and using equation (B36), or by 
the amplitude ratio and using equation (B37). 

The convective heat -transfer parameters calculated as a function of phase 
lag <p, wall-properties-thickness -frequency parameter rjli, wall temperature location 
x/L, and heat-transfer coefficient ratio R, using equation (B36), is plotted in figure 2. 

A second plot of the convective heat -transfer parameters calculated as a function of 
wall temperature amplitude ratio wall-properties -thickness -frequency para- 

meter 7]li, wall temperature location x/L, and heat -transfer coefficient ratio R, using 
equation (B37), is shown in figure 3. The R = 0 charts are included because they show 
the special case of an insulated surface at x/L = 1.0, that is, at h^ = 0 (ref. 4). In 
figure 2 the charts for x/L =1.0 have curves of \p against cp which appear to give 
negative values for p. This is explained by letting h^ approach infinity so that 
also approaches infinity. As this happens, R approaches zero, which forces <p to 
approach zero provii^ that positive values of p can exist for p = 0. Figure 3 shows 
that for x/L =1.0 the curves for p are double valued. The choice of values can be 
made using the phase lag angle solution of figure 2 to pick the correct value. 

Figures 2 and 3 are included in this report to aid the experimentalist in (1) determin- 
ing the frequency and amplitude of the hot -gas temperature oscillation that is required to 
give a measurable wall temperature response, and (2) calculating the convective heat- 
transfer coefficient. When using the charts to determine wall temperature response, a 
calculation should be made to determine whether the coolant temperature will oscillate 
significantly, which would violate the fourth assumption used in the derivation of the solu- 
tion. An example is given here to show how the heat-transfer coefficients may be calcu- 
lated using figure 2. 

(1) Assume for the moment that the temperature drop across the wall is small. Then 
the measured wall temperature can be taken equal to T(0) and T(L). The mean value of 
the hot fluid Tq, wall T(0), and coolant T^ temperatures are substituted in equation (9) 
and the heat -transfer coefficient ratio R is calculated. The value of R and the location 
of the measured wall temperature, x/L =0 or 1. 0 for these charts, determines which 
chart in figure 2 is used. 
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Convective heat -transfer parameter, = (Htj/Kq) 
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0 10 20 30 40 50 60 70 80 


Phase lag angle, cp , deg 

(a) Wall temperature measurement location, x/L = 0; ratio of convective 
heat-transfer coefficients, 0 (insulated wall) . 

Figure 2. - Convective heat -transfer parameter as function of phase lag 
angle. 


10 



Convective heat -transfer parameter, i//q = {Kri/hQ^ 





Convective heat -transfer parameter, xj^Q = (Kti/Kq) 
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0 10 20 30 40 50 60 70 80 

Phase lag angle, w, deg 


(c) Wall temperature measurement location, x/L = 0; ratio of convective 
heat -transfer coefficients, 0.5. 

Figure 2. - Continued. 
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Convective heat -transfer parameter, = (K 7 ?/hQ) 
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Phase lag angle, cp, deg 


(e) Wall temperature measurement location, x/L = 0; ratio of convective heat 
transfer coefficients, 1.0 


Figure 2. - Continued. 
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Phase lag angle, <p, deg 


(i) Wall temperature measurement location, x/L = 0; ratio of con 
vective heat -transfer coefficients, 3.0. 


Figure 2. - Continued. 
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( 3 ) Wall temperature measurement location, x/L = 1.0; ratio of convective heat- 
transfer coefficients, 3.0. 


Figure 2. - Continued. 
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0 10 20 30 40 50 60 

Phase lag angle, <p, deg 


(k) Wall temperature measurement location, x/L = 0; 
ratio of convective heat -transfer coefficients, 5.0. 

Figure 2. - Continued. 
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Convective heat -transfer parameter, = (KTi/hg) 



Phase lag angle, (p, deg 

(1) Wall temperature measurement location, x/L = 1.0; ratio of convective heat- 
transfer coefficients, 5.0. 

Figure 2. - Continued. 
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Convective heat -transfer parameter, = (Krj/h.„) 



0 10 20 30 40 50 

Phase lag angle, cp, deg 


(m) Wall temperature measurement location, x/L = 0; 
ratio of convective heat -transfer coefficients, 7.0. 


Figure 2. - Continued. 
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0 10 20 30 40 50 60 


Phase lag angle, <p, deg 

(o) Wall temperature measurement location, x/l. = 0; ra- 
tio of convective heat -transfer coefficients, 10.0. 

Figure 2. - Continued. 
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0 10 20 30 40 50 60 70 80 90 

Phase lag angle, cp, deg 


(p) Wall temperature measurement location, x/L = 1.0; ratio of convective heat- 
transfer coefficients, 10.0. 


Figure 2. - Concluded. 
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(a) Wall temperature measurement location, x/L - 0; ratio of convective heat -transfer coei- 
ficient (insulated wall), 0. 


Figure 3. - Convective heat-transfer parameter as function of wall temperature amplitude 
rat’O. 
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Convective heat -transfer parameter, = (K? 7 /hp) 
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(b) Wall temperature measurement location, x/L = 1. 0; ratio of convective heat -transfer 
coefficients, 0 (insulated wall) . 

Figure 3. - Continued. 
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(c) Wall temperature measurement location, x/L = 0; ratio of convective heat -transfer co- 
efficients, 0. 5. 


Figure 3. - Continued. 
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Convective heat -transfer parameter, \j/Q = (Krj/hfj^ 
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Amplitude ratio, 

(d) Wall temperature measurement location, x/L = 1.0; ratio of 
convective heat -transfer coefficients, 0. 5. 

Figure 3. - Continued. 
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Convective heat -transfer paramter, = {Krj/h^ 
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(e) Wall temperature measurement location, x/L = 0; ratio of convective heat -transfer coef- 
ficients, 1.0. 

Figure 3. - Continued. 
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Convective heat -transfer parameter, = {Kr]/h^ 
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(f) Wall temperature measurement location, x/L = 1.0; 
ratio of convective heat -transfer coefficients, 1.0. 


Figure 3. - Continued. 
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(h) Wall temperature measurement location, x/L = 1.0; 
ratio of convective heat -transfer coefficients, 2.0. 


Figure 3. - Continued 
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Convective heat -transfer parameter, = (Kt7Aq) 




Convective heat -transfer parameter, \j/„ = {Kr]/h„) 
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Amplitude ratio, 

(k) Wall temperature measurement location, x/L = 0; ratio of convective heat -transfer co 
efficients, 5.0. 


Figure 3. - Continued. 



Amplitude ratio, 


(1) Wall temperature measurement location, x/L = 1.0; ratio of convective heat -transfer co- 
efficients, 5.0. 

Figure 3. - Continued. 
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(m) Wall temperature 
efficients, 7.0. 


Amplitude ratio, 

sasurement location, x/L = 0; ratio of convective heat -transfer co 


Figure 3. - Continued. 
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(n) Wall temperature measurement location, x/L = 1.0; ratio of convective heat -transfer co- 
efficients, 7.0. 


Figure 3. - Continued. 
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(o) Wall temperature measurement location, x/L = 0; ratio of convective heat -transfer co- 
efficients, 10.0. 


Figure 3. - Continued. 


40 




Amplitude ratio, 


(p) Wall temperature measurement location, x/L = 1.0; ratio of convective heat -transfer 
coefficients, 10.0. 


Figure 3. - Concluded. 
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(2) The frequency of the temperature oscillation, density, specific heat, thermal con- 
ductivity, and wall thickness are used to calculate the value of tjIj, This determines which 
curve is to be used on the selected chart. 

(3) The phase iag (in deg) of the wall temperature behind the forced fluid temperature 
is determined experimentally from simultaneous oscillograph traces of these tempera- 
tures. Entering the selected chart with the phase lag angle on the appropriate rjL curve 
gives a value for the convective heat -transfer parameter i//q. 

(4) The convective heat -transfer coefficient on the hot-gas side is calculated by sub- 
stituting ij/Q from step (3), the frequency of the temperature oscillation, and the wall 
properties into equation (7). 

(5) The convective heat -transfer coefficient on the coolant side is calculated using 
either equation (8) or more simply from: h^ = Rh^. 

Obtaining the correct value of R in figures 2 and 3 may present some difficulty, 
since only one wall temperature is assumed to be measured and since equation (9) calls 
for both T(0) and T(L). This difficulty, however, may be circumvented as explained as 
follows. 

If the temperature differences in equation (9) are large enough that the temperature 
drop across the wall does not significantly affect the ratio R, the required relation 
between and is given by equations (8) and (9) by assuming that T(0) = T(L) 

K R is affected significantly in equation (9) by the wall temperatures T(0) and 
T(L) not being equal, an iternative process may be used. The initial value for R is 
calculated from equation (9) assumii^ that T(0) equals T(L). Resulting values for h^ 
and h^ will allow a heat flux q to be calculated using the measured wall temperature 
(assumed for this illustration to be at x = 0, although any known location in the wall will 
also be soluble) and the following equation; 

q = h^(T^ - T(0)) (10) 

The temperature drop through the wall is given by 

T(0) - T(L) = ^ (11) 

K 


so that 


T(L) = T(0) - ^ 
K 


(12) 


Use of the new value of 


T(L) in equation (9) gives an improved value for 


R, which. 
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in turn, can be used in the transient solution to obtain improved values of heat-transfer 
coefficients. The process is repeated until the heat -transfer coefficients no longer 
change significantly. It has been found that between three to six interations will match 
the heat convected to the wall to that which is conducted throi^h the wall and convected 
from the wall. The agreement is within 1 percent. 

Since charts cannot fit exactly the conditions of an experiment without interpolations, 
equation (B26) has been rearranged as shown in appendix B (eq. (B51)) to give a more 
direct solution for This equation involves the convective heat -transfer coefficient 

ratio R, along with rjL, x/L and the phase lag angle cp, and it can be solved using the 
classic cubic equation method. The correct choice of root for equation (B51) is the 
positive, real root. 

Equation (B27) has not been rearranged to give or in terms of the ampli- 
tude ratio, since the amplitude ratio requires higher frequencies to get measurable am- 
plitude ratio changes than the phase lag equation. Also, using the phase lag eliminates 
the need for having the absolute value of ATq and thereby permits the use of a trace of 
chamber pressure in the case of a rocket engine to replace the oscillating trace of the 
hot -fluid temperature. Use of the phase lag equation also eliminates the double value 
problem encountered with the amplitude ratio equations. However, this equation has been 
programmed and gives h^ and h^ as functions of amplitude ratio. Both equations 
(B51) and (B27) have been programmed (unpublished report by G. Aling and R. Huff) in 
FORTRAN IV and are operational on the IBM 7094-2/7044 direct-coupled system. The 
titles of the programs are JCL and BL, respectively. 


Lumped Wall Properties (Case 2) 

The assumptions used in the case 1 solution have been applied to the case 2 solutions 
with one exception. Assumption (1) has been eliminated, and in its place it has been 
assumed that the temperature gradient through the wall is small or that the wall tempera- 
ture is measured at a point in the wall where the thermal properties are considered to 
be lumped (appendix C). 

The governing differential equation can be written by equating the difference between 
the heat transferred to the wall and the heat that leaves the wall to the heat stored in the 
wall. In writing the differential equation which describes this condition, it is possible to 
account for the thermal conductivity of the wall by using a modified overall heat -transfer 
coefficient. The form of the differential equation is the same as that which results when 
the thermal resistance of the wall is neglected. Therefore, the differential equation that 
neglects the effects of thermal conductivity is used in the initial derivation and then the 
solution is modified to account for the thermal conductivity of the wall material. The 
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differential equation used in the initial derivation (appendix C) is 


pcL = h^{T^ - T) - h^(T - T^) (13) 

9t 

The hot -gas -side fluid temperature Tg is ariven sinusoidally and is given by 

Tq = sin o)T (14) 

The solution to equation (13), using Laplace transform techniques is detailed in 
appendix C and is given here in functional form as 


Ti(t) 


__ h„ + h 

"^c’ ^'^G’ P’ ^G’ ^ 




pcL 


Starting transient 




N(Tq, T^, h^_, hj 




Offset from conditions before time zero 


^i(ATg, CO, p, c, L, hj^, hp)sin(coT - <Pj) 
Steady -state oscillation 


where 


cp^= (pip, c, L, CO, h^, h^) 

Laplace transforms are used because they yield a solution that gives the wall temperature 
history from just before the initial gas temperature fluctuation to the time when the wall 
temperature reaches a steady -state (neutral) oscillation. The solution for the wall tem- 
perature is given by equation (C9). The phase lag is given by equation (CIO). 

By assuming that the ratio of the convective heat -transfer coefficients R can be 
calculated from equation (9), the equations for the heat -transfer coefficients in terms of 
either phase lag angle or amplitude ratio can be derived from equations (CIO) and (Cll), 
respectively. These are 



pcLct; 

(1 + R)tan cp^ 


( 15 ) 



and 


- 1/2 


- pcLo) 


(1 + R)^ 


The coolant side coefficient from the definition of R is 


( 16 ) 


^ ( 17 ) 

The time required for the decay of the initial transient of the wall temperature so that 
it is influenced only by the hot -gas temperature oscillations can be determined by using 
the first three terms of equation (C9) (starting transient). This calculation serves only 
as an approximation as to the time required to reach a steady -state oscillatii^ condition. 
The approximation should be good when the thermal dtffusivity of the wall material is 
high. 

To account for the thermal conductivity of the wall material when using the lumped- 
wall -property method, equation (13) can be modified by replacing the convective heat- 
transfer coefficients with a modified overall coefficient of heat transfer. The details of 
this procedure are given in the second section of appendix C. This modification results 
in the following equations for the wall temperature amplitude ratio and phase lag angle 

From equation (CIS) the amplitude ratio is 



From equation (C23) the phase lag angle is 


^2 “ 


tan 





L xpj 


r/w!G3l\ 
^ , \ 

IH- 

V L/ 


(19) 
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Equations (18) and (19) can be rearranged so that the convective heat -transfer para- 
meter can be calculated directly. Equation (18) becomes 


77LR 




G2 


■(-?) 


cos 


(P2 


^m2 


AT 


GJ 


‘^m2 

AT^ 


(1 + R) - cos ^2 


( 20 ) 


and equation (19) becomes 




G2 


= il±^tan 

4?7L 


cp2 


rfL 

2 


^+11 


(1 + H) 

27jL 


R 



tan ( 


CP 2 } -4 


J 




A 1/2 


II - 


X 


R - 5. tan <p^ 
L / 2 


(21) 


COMPARISON OF DISTRIBUTED (CASE 1) AND LUMPED (CASE 2) SOLUTIONS 

Generally, case 2 solution not accounting for conductivity is preferred because of its 
mathematical simplicity; however, large errors in the convective heat -transfer coeffi- 
cients may be introduced by using the case 2 solution. 

The best way to determine whether the case 2 solution can be used in a particular 
experiment is to compare the heat -transfer coefficients calculated using both the case 1 
and case 2 solutions over the expected range of operating conditions. If the difference is 
negligible, the case 2 solution would, for simplicity purposes, be used. 

To compare case 1 with case 2, a plot of the revised amplitude ratio (fig. 4(a)) and 
phase lag angle (fig. 4(b)) versus the tangent of the case 2 phase lag angle is shown in 
figure 4. The case 2 solution provides a base line, for comparison purposes on both 
plots, that does not shift with changing tjL, x/L, or R values. The reason for the non- 
shifting base line in figure 4(a) is that the revised amplitude ratio is defined as 
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Revised amplitude ratio, 



Tangent of case 2 phase lag angle accounting for wall thermal conductivity, tan ‘P 2 
(a) Wall temperature amplitude ratio. 

Figure 4. - Comparison of case 1 and case 2 solutions; ratio of convective heat-transfer coefficients, 5.0. 
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accounting for wall thermal conductivity, tan f g 


Phase lag angle. 


4. - Concluded. 



which, from equation (18), is equal to cos Plotting the cos <P 2 versus the tan <p 2 
gives the base Line curve shown in figure 4(a), which is not a function of r]L, x/L, or R. 
The reason for the nonshifting base line in figure 4(b) is, from inspection of equation (19), 
that the base line is merely a plot of the phase lag angle versus the tangent of the phase 
lag angle. The amplitude ratio (fig. 4(a)) and the phase lag angle (fig. 4(b)) are presented 
in this manner so the reader may relate these solutions to their electromechanical analo- 
gies, which are often plotted this way. The base lines are the same as those which re- 
sult from analysis of an electrical series resistance -capacitance circuit driven by a 
sinusoidal voltage. 

In order to compare the amplitude ratio for the case 1 solution with the base line, the 
case 1 solution for amplitude ratio 0 ^/ATq (eq. (B27)) was multiplied by the bracketed 
quantity in equation (22), and this revised plotted on figure 4(a) versus the 

tan q >2 calculated using equation (19). In this way, the case 2 solution accounting for 
thermal conductivity can be used when the case 1 plots fall on or near the base line curves 
of case 2. 

Similarly, the phase lag angles calculated using case 1 (eq. (B26)) may be compared 
in figure 4(b) to the base line of case 2. The case 1 phase lag angles are plotted versus 
the tan <P 2 , calculated using equation (19). Again, when the curves fall close to or on the 
base line the case 2 solution accounting for the thermal conductivity can be used in place 
of case 1. 

The set of curves shown in figure 4 are calculated for R = 5. 0 and x/L = 0 and 1. 0. 
These conditions only affect case 1 curves, because the base-line curves are not influenced 
by these parameters. The general shape of the curves in figure 4 is the same for values 
of R between 0. 5 and 10. 0. Inspection of figures 4(a) and (b) shows that small values of 
77L tend to allow the use of case 2 solutions. 

To show the order of magnitude of the errors in convective heat-transfer coefficients 
that could result by using the case 2 solution instead of case 1, two examples will be pre- 
sented (based on phase lag angle calculations). The example shown as a function of fre- 
quency in figure 5 assumes a rocket engine having a stainless -steel wall 0. 012 inch 
(0. 3048 mm) thick and a hot-gas-side heat-transfer coefficient of 0. 0022 Btu per square 
inch per °R (6.472 kW/(m^)(K)). 

For the rocket engine example, an R value of 6. 1 is considered representative; 
however, a range of R values are presented. For these conditions and an x/L =1.0 
(fig. 5(a)), errors from 30 percent at 6 hertz {r]h = 0. 606) to 20 percent at 1. 0 hertz 
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(a) Wall temperature maasurement location, x/L =1.0 and 0. 5. 

Figure 5. - Percent error in convective heat-transfer coefficients using 
lumped wall properties for rocket engine example. Assumed wall 

q 

conditions: density, 0. 286 pound mass per cubic inch (7. 92 g/m ); 
specific heat, 0. 1358 Btu per pound mass per °R (1. 5686x10^ J/(kg)(K)); 
thermal conductivity, 0. 287x10 ^ Btu per inch per second per °R 
(21. 46 J/(m)(sec)(K)); thickness, 0. 012 inch (0. 3048x10"^ m); gas-side 
heat-transfer coefficient, 0. 0022 Btu per square inch per second per °R 
(6.472 kW/(m^)(K)); h^ = R h^. 

(tjL = 0. 25) are shown to be possible. At x/L = 0 (fig. 5(b)), errors are reduced to 
-22 percent at 6 hertz and -2 percent at 1. 0 hertz. 

The percent errors are based on the assumed value of h. The case 1 solution 
(eq. (B26)) was used to calculate the phase lag angle; this angle was then used in the 
case 2 solution (eq. (15), neglecting the thermal conductivity of the wall) to calculate the 
heat -transfer coefficient. 

The second example assumes a much lower heat-transfer coefficient on the hot -gas 
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(b) Wall temperature measurement location, x/L = 0. 


Figure 5. - Concluded. 


side (by a factor of 10), considered to be representative of a heat exchanger. Figure 6 
shows coefficient errors for the heat exchanger example. This example assumes a 
stainless -steel wall 0.06 inch (1. 524 mm) thick, and a hot-gas-side heat-transfer coeffi- 
cient of 0. 00021 Btu per square inch per second per °R (0. 61782 kW/(m^)(K)). An R 
value of 1. 16 was considered representative of the heat ex charter case; again however, 
a range of R values are presented. The coefficient errors were calculated using the 
same procedure followed for the rocket engine example. Under these conditions for 
x/L = 1. 0 (fig. 6(a)), the errors were as large as 62 percent at 0. 2 hertz (tjL = 0. 566) 
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Error in using case 2 to calculate convective heat-transfer 
coefficient, (h - h „)/h x 100, percent 


x/L 


1.0 

I Frequency at which phase 
lag reaches 20° 


Ratio of convective 
heat -transfer 
coefficient, 



(a) Wall temperature measurment location, x/L =1.0 and 0. 5. 

Figure 6. - Percent error in convective heat -transfer coefficients using 
lumped wall properties for heat exchanger example. ■ Assumed wall con- 
ditions: density, 0.286 pound mass per cubic inch (7.92 g/m ); specific 
heat, 0. 120 Btu per pound mass per °R (0. 5024x10^ j/kg)(K)); thermal 
conductivity, 0.2411x10^ Btu per inch per second per °R (18.02 
j/(m)(sec)(K)); thickness, 0.06 inch (1.524 mm); hot-gas side convec- 
tive heat -transfer coefficient, 0.00021 Btu per square inch per second 
per °R (0. 6178 kW/(m^)(K)); h^ = R h^^. 
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Error in using case 2 to calculate convective heat-transfer 
coefficient, (h - o)A ^ percent 

CaSc Ci 



Figure 6. - Concluded. 
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and as small as 3 percent of 0. 01 hertz (t]L = 0. 1267). For x/L = 0 (fig. 6(b)), the errors 
exceed -100 percent at 0. 2 hertz (?7L = 0. 0566) but are as small as -3 percent at 0.01 hertz 
(r?L = 0. 1267). 

Both examples show that decreasing the frequency (also tjL) causes the errors to 
decrease. However, in the rocket engine example for x/L = 1. 0 the errors remained at 
20 percent, showing tiiat low values of frequency or t/L will not necessarily give negli- 
gible differences between case 1 and case 2 solutions. This is a result of the difference 
between the assumptions used in formulating the case 1 and case 2 solutions. At low 
values of tjL, R values near 1. 0 yield lowest errors. 

To compare the two solutions accounting for thermal conductivity of the wall material, 
equation (21) was used to calculate the convective heat-transfer parameter with the phase 
lag angle obtained from the case 1 solution (eq. (B26)). The results for the rocket engines 
example but for an R = 1. 0, are shown in figure 5 by the dashed lines. For the rocket 
engine conditions, the heat -transfer coefficients obtained from equation (21) (considering 
conductivity) instead of equation (15) are not significantly improved. 

Differences between the case 1 and case 2 solutions do exist even if the frequency is 
allowed to approach zero. This has been shown in figures 5 and 6 and is expected because 
of the differences in the assumptions used to generate the two solutions. Generally, 
though it appears that a combination of a low value of rjL and R is required to minimize 
the differences between the heat-transfer coefficients calculated using the case 1 and 
case 2 solutions involving phase lag angles. 


CASE 2 SOLUTION WITH MIDPOINT WALL TEMPERATURE MEASUREMENT 

An attractive possibility for reducing the differences between case 1 and case 2 
solutions is to measure the wail temperature at the center of the v/all, x/L = 0. 5. This 
will obtain a phase lag angle which should be nearly equal to the average phase lag angle. 
This results because the phase lag angle is approximately proportional to the distance 
through the wall tliat the wave has traveled. In figure 4(b), it can be seen that the averages 
of x/L = 1. 0 and x/L = 0 for the same rjL values are always much closer to the 
case 2 base line than either surface alone is. If the v/all temperature is measured at the 
center, the case 2 solution, in terms of phase lag angle, can possibly be used to calculate 
the coefficients. 

The curves given in figure 5(a) for x/L = 0. 5 show that under certain conditions the 
errors caused by using the case 2 solution to calculate the coefficients tend to be reduced. 
This is true when the comparison is made between the x/L =1.0 and 0. 5 case. How- 
ever, a comparison of figure 5(b), x/L = 0, R = 10 to the x/L = 0. 5 curves of figure 5(a) 
shows that the errors in coefficient are not reduced. Therefore, it can only be stated that. 
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when the wall temperature is measured at the midpoint, the errors due to using the case 2 
solution to calculate the coefficients tend to be reduced. 


OPTIMUM PHASE LAG ANGLE 

For the insulated wall case reference 1, an optimum range of phase lag angles exists. 
The percent error in the coefficient due to an error in measuring the phase lag angle is 
minimized in this range. The same results are obtained when using the case 2 solution 
given in this report. The equation for an incremental change in convective heat -transfer 
coefficient is obtained from equation (15) and is 

? (23) 

h dcp sin 2(p 

As shown in reference 2, the minimum value exists a.t cp = 45° (0. 785 rad). Further- 
more, cp can range from 20° to 80° (0. 349 to 1. 396 rad) before the errors in h begin to 
increase greatly. Hence, in using the simpler case 2 solution, a guide line is established 
which shows that a minimum phase lag angle of 20° (0.349 rad) and a maximum of 80° 

(1.396 rad) are desirable (to minimize errors in h due to inaccuracies in measuring the 
phase lag angle). However, it appears from inspection of figure 2 that the case 1 solution 
imposes an upper limit on the phase lag which is dependent on both x/L and R. The upper 
limit is best determined by inspection of figure 2 by examining the slope of the i// versus 
cp curves. When the Ax{//Acp becomes large enough so that the error in the measurement 
cp gives an unacceptable error in xj/ a lower frequency must be used. This will give 
smaller phase lags and better accuracy. 

In figures 5 and 6 the frequencies at which a phase lag angle of 20° is reached are 
indicated by arrows shown on each of the curves. These indicate the start of the optimum 
frequency range. Higher frequencies will increase the phase lag angle. As can be seen, 
these frequencies are low enough (below 3 Hz) so that they should be easily obtainable with 
current control systems. 


SOLUTION OF PROBLEM WHEN COOLANT TEMPERATURE 
IS DRIVEN SINUSOIDALLY 

The preceding solutions have been derived for the case in which the hot -gas tempera- 
ture is driven sinusoidally. The extension of this technique to the case in which the 
coolant temperature is driven requires the followir^ substitutions to be made: 
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he 

R 


T = T ’ 


G 


Ill III 

'^case 1 

= h0 

r C 1 

= ^gJ 

|.case 2 

h' 

c 

T(0) - T^ 

^G 

Tq - T(L) 

T = ' 

K - 


For the case 1 and case 2 solutions these substitutions change the stated location of the 
coolant -wall surface to x = 0 and that of the hot-gas wall surface to x = L (i. e. , x = 0 
at the surface on which the fluid temperature is oscillated). Making these substitutions 
allows the charts in figures 2 and 3 to be used to calculate the heat-transfer coefficients 
when the coolant temperature is varied sinusoidally instead of the hot -gas temperature. 
When making these substitutions in figures 2 and 3 the negative sign of the coefficient is 
discarded and only the absolute value of the convective heat -transfer parameter \p is 
used. 


SUMMARY OF RESULTS 

An analytical investigation has been performed in which the temperature response of 
a convectively heated and cooled plate to a sinusoidally varying driving temperature is 
determined. The temperature response of the wall can be used to determine both con- 
vective heat -transfer coefficients. Two independent methods for calculatir^ the coeffi- 
cients are presented. One uses the phase angle of the measured wall temperature lag 
behind the driven (hot gas or coolant) temperature. The other method uses the ratio of 
the amplitude of the wall temperature to the amplitude of the driven temperature. 

Each of these methods requires two independent relations between the coefficients 
in order to solve for the coefficients. The first relation is knowledge of the sinusoidal 
variation of the driven fluid and wall temperatures with time, so that the phase lag angle 
between these quantities (or amplitude ratio) can be determined. The second requires the 
time averaged values of the driven, wall, and coolant temperatures. These are used to 
calculate the convective heat-transfer coefficient ratio. These two relations provide the 
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information needed to calculate the coefficients. 

Two solutions for the wall temperature response are presented. The first, case 1, 
uses distributed wall properties and accounts for low thermal diffusivities and thick walls. 
In this case, the location of the measured wall temperature is important and the solutions 
for the convective heat -transfer coefficients in terms of phase lag angle or amplitude 
ratio are lengthy. Charts giving the solutions for the coefficients over a selected range 
of parameters are included. The other solution, case 2, assumes the wall has a lumped 
capacitance and thermal conductivity. In this case, the solutions for the convective heat- 
transfer coefficients using either phase lag angle or amplitude ratio are simple. The 
solutions are analyzed by comparing the convective heat -transfer coefficients calculated 
using the distributed solution to those calculated using the lumped solution. 

Measuring the wall temperature at the midpoint in the wall tends to minimize the 
errors in heat-transfer coefficients when the lumped wall property solution is used to 
calculate them. 

An optimum range of phase lag angles exists within which minimum errors in con- 
vective heat-transfer coefficients result from errors in measurement of phase lag angles. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 5, 1969, 

122-29. 


57 


APPENDIX A 


SYMBOLS 


c specific heat of wall material, 

Btu/(lbm)(°R); J/(kg)(K) 

e natural logarithm base, 

2. 71828 

F(t) part of product solution for 

0(x, t) which is independent 
of X 


f frequency of temperature os- 

cillation; Hz 

h time averaged convective 

heat -transfer coefficient in- 
dependent of the impressed 
sinusoidal temperature vari- 
ation, Btu/(in.^)(sec)(°R); 
W/(m2)(K) 

i imaginary number, yTi 

K thermal conductivity of wall 

material, Btu/(in.)(sec)(°R); 
J/(m)(sec)(K) 

L thickness of wall material, 

in, ; m 

jg’[T] Laplace transform of tem- 

perature 

M, N, P, V functional notation 

q heat flux per unit area, 

Btu/(in.^)(sec); W/(m^) 


R ratio of convective heat- 

transfer coefficients, h^/h^ 

S Laplace transform variable 


T 


AT 

T 


T(L, t) 
T(0, t) 
T(x, t) 


Ts(x) 

U 


u(t) 

X(x) 


X 


temperature (temperature of 
wall when not subscripted), 
°R; K 

half amplitude of temperature 
oscillation, °R; K 

average or mean value of tem- 
perature, °R; K 

wall temperature, function of 
time at x = L, °R; K 

wall temperature, function of 
time at x = 0, *^R; K 

wall temperature, function of 
distance (measured into 
wall from hot surface) and 
time, °R; K 

steady -state part of T(x, t), 
°R; K 

modified overall heat -transfer 
coefficient which accounts 
for the thermal conductivity 
of the wall material, 

Btu/(in. ^) (se c) (°R) ; 
W/(m^)(K) 

unit step function 

part of the product solution 
for 0(x, t) which is inde- 
pendent of time, °R; K 

distance measured from the 
heated surface of the wall 
into wall toward cooled 
surface, in. ; m 
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Oi 

V 

e 

p 

T 

(P 

(jO 

Subscripts: 

c 

cw 

G 

GW 


thermal diffusivity, K/pc, 
in.^/sec; m^/sec 

frequency and wall property 
parameter, '^u>/2a, 1/in. ; 
1/m 

transient part of appropriate 
wall temperature solution, 
°R; K 

density of wall material, 
Ibm/in.^; g/m^ 

time, sec 

phase lag angle between the 
driving gas temperature 
and the responding wall 
temperature, deg; rad 

convective heat -transfer 
parameter Kri/h, dimen- 
sionless 

angular velocity of tempera- 
ture oscillation, 27rf rad/sec 


coolant 

coolant side of wall, x = L 

hot gas or fluid 

hot -gas side of wall, x = 0 


m maximum value 

r revised by using 

first-order solution to cor- 
rect the value to range 
from 0 to 1. 0. 

case 2 refers to values calculated 

usit^ the appropriate 
case 2 solutions 

0 refers to conditions before 

the start of the tempera- 
ture oscillation 

1 case 2 first-order differen- 

tial equation solution not 
accounting for the thermal 
conductivity and using 
lumped system 

2 case 2 first-order differen- 

tial equation solution 
accounting for the thermal 
conductivity using modified 
overall heat -transfer 
coefficient 

Superscript: 

' solutions in which the coolant 

temperature is driven 
sinusoidally. In appendix B 
the prime indicates a deri- 
vative 
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APPENDIX B 


DERIVATION OF THE CASE 1 SOLUTIONS - DISTRIBUTED WALL PROPERTIES 

Figure 1(a) shows a wall through which heat is being transferred from one fluid to 
another. The heat is transferred across the fluid-wall boundaries convectively. It is 
required that the convective heat-transfer coefficients at the plate surfaces be deter- 
mined experimentally by measuring the two fluid temperatures and only one wall temper- 
ature. This can be accomplished with the given information by oscillating one of the fluid 
temperatures sinusoidally and measuring either the phase lag angle between the fluid and 
wall temperature or the amplitudes of the fluid and wall temperature. A sketch showing 
what the wall temperature might look like is given in figure 1(b) . The solutions for the 
coefficients are found by first finding the equation for the wall temperature response to 
the sinusoidally driven fluid temperature and then solving this equation in terms of the 
known quantities for the coefficients. 


General Solution 

The response of a convectively heated and cooled wall to a sinusoidally varying hot- 
gas temperature is derived here. The solution accounts for the thermal properties of 
the wall material being distributed throughout the wall. The following assumptions are 
used in the derivation: 

(1) The heat flows through the plate in the x-direction only (one -dimensional heat 
conduction) . 

(2) The wall or plate properties (density, specific heat, and thermal conductivity) 
are constant. 

(3) The convective heat -transfer coefficients on both sides of the wall are constant. 

(4) The coolant temperature is constant. 

The second -order partial differential equation for transient, one -dimensional, heat 
conduction in a plate is solved. The boundary conditions are that the plate is convectively 
heated on one side by a fluid having a sinusoidally varying temperature and cooled on 
the other side by a fluid which is not significantly affected by the oscillating heat flux 
from the wall. The solution follows. 

The differential equation is 

^=1”; (Bi) 

ax^ « 


60 



where a = K/pc, the thermal diffusivity. 
The boundary conditions are 


ho [T(j - T(0, r)] 

h, [t(L, t) - T J 

The driving temperature is 



at X = 0 

(B2) 

dx 


_ 3T(L, t) 

at X = L 

(B3) 


8x 


Tq = AT^ e 


-IWT . 7 ^ 


+ T, 


G 


(B4) 


Assume the solution has the form of a transient plus a steady -state value then 

e(x, t) = T(x, t) - Tg(x) (B5) 

where Tg is the steady-state temperature distribution prior to driving Tq(t) and T(x,t) 
is the temperature response of the wall to the harmonically driven hot -gas temperature. 
The solution for g{x, t) follows. 

Using boundary condition equation (B2) and substituting for Tq the driving gas tempera- 
ture, equation (B4), yields, after subtracting the steady -state heat transfer 


ATq e - 0(0, t) 


K 30(0, r) 
hQ 3x 


Boundary condition equation (B3) can be rewritten as 
The differential equation becomes 

3^0(x, t) _ ^ 30(x, t) 



(B6) 


(B7) 


(B8) 


Assume a product solution then 
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0 = X(x) • F(r) 


(B9) 


F(t)X’(x) 

0X 


X(x)F'(t) 

9t 

2 

F(r)X”(x) 

2 

3x 

Substituting this into differential equation (B8) gives 

F(t)X”(x) = 1x(x)F(t) 
a 

X»'(x) _ 1 F'(r) 
X(x) a F(t) 


(BIO) 


Since either side of equation (B9) is independent of the variable appearing on the other 
side, it is concluded that either side must be equal to a constant. This constant is 

p 

assumed to be X and from equation (B9) the following equations have been written: 

X" - X^ (Bll) 


F’ = X^cvF (B12) 

Three choices for X^ are possible: X^ <0, X^ = 0, X^ > 0. If X^ > 0, the solution will 

2 2 

not be periodic in x: therefore, X > 0 is rejected. If X = 0, the solution will not be 

2 2 

periodic; therefore X = 0 is rejected. If X < 0, a solution exist and is found as 
follows: 

From equation (Bll) 


X” + X^ = 0 

The characteristic is 

m^ + X^ = 0 
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From which 


m = iX 


From Wiley (ref. 5, p. 88) 


X = Cj e + <^2 6 

p 

From the derivation for F given below in which the equation for X is deduced, 
lowing equations for X are derived; 


a 


X = 


X = ± 



X = ± (1 + i) 



Then 


X = c X ^ gT(l+i) Vgo/2q'x 


From equation (B12) 


F 


a 


In F = 


- X ar 


F = 


-X" 


c e 


Q!T 


the fol- 


(B13) 
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2 — 

If A is real, F cannot be periodic since 

2 

g-X ar _ gQgjj x^OT - sinh X^ar 

2 

But if X is imaginary, then 


e"iX = cos X^Q!T - i sin X^ot 
2 

which is periodic as required. If X = io)/ a, where cu is the angular velocity. 
2 

X ar terms in this equation reduce to the usual wt form. Hence 


F = 


c e 


-iWT 


Substituting equations (B13) and (B14) into (B9) gives 


0 = c e 


-10)T 


cie 


±i V(a)/2a)xT'\/(a)/2cv)x Ti V(o)/2a)x± V(cu/2o’)x 

+ ^2 


The choice of signs in the exponent is arbitrary. Combining the constants (cc^ = 
CC2 = B) and rewriting this equation gives 




where 7j = Va;/2a. 

Application of the boundary conditions requires the 90/9x. This is given by 

^ = A ^ ^ B g77X-i(cuT+rjx) 

9x 

= 77 A(i - i)e-^x-i(wT-7,x) ^ ^ B(1 - De^^-H^r+rjx) 

9x 


The 

(B14) 

(B15) 
A and 

(B 16 ) 

(B17) 

(B18) 
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Now equation (B6) gives 


AT^ -0(O,r)=^MO^ 

ax 

AT^ e-i^’- - A e-i^^ - B = - ?^ (A - B)(i - 

"g 

AT„ - A- B = ?^(B- A)(i - 1) 

° >’G 



- A 

(i - 1) - 1 

'’g 


_^G 


Applying the boundary condition equation (B7) at x = L gives 


(B19) 


0(L,r)=-KMkj) 


0x 


A e e wt-tjL) g ^ -i( wt+?]L) ^ ^ ^ -tjL -i( wt -77L) _ ^ ^ ^t/L - i( wr+rjL) 


A 


1 + M (i - 1) 

g-?7L-i(wT-7]L) ^ _g 

1 

+ 

_ 


L h, J 


^77L-i(o)T+77L) 


A 

B 


1 + ^^(1 -i) 

he 


^27)L-i27jL 


(B20) 


1 + (i - 1) 
h„ 


Using this expression and substitutii^ in equation (B19), which was obtained from bound- 
ary condition equation (B6), the equations for the constants A and B have been ob- 
tained as follows: 
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AT =BM(i-l) + l +B 




1 + M(i_i) 


AT^ 1 + 5!I (i 

L '■c 


(i - 1) - 1 e 


^2r]L-i2r]L 


i - 1) 


^(i-l) + l l + + I^(i_i) _i e 2 ’ 7 L-i 2 T?L (^ 2 ^^ 


-AT^ 1 + M (1 - i) e2r]L-i2?7L 

°L •'c I 


52 (i - 1) + l1 + !^ (i - 1)1 + fl + 52 (1 - M (1 - 1) - x1e22L-12ijL (222) 


Working with the common denominator in the equations for the constants A and B (eqs. 
(B21) and (B22)) and defining x}^ =Kri/h results in 

- 1 ) + ^ |Vj.(i - 1 ) + ^ + - 1 ) + ^ [VG(i - 1 ) - ^ g 2 rjL-i 2 ? 7 L 

|jl - \Pq) + ii//J [(I - xp^) + + Iji + - ip^ [1 (1 + ;^^) + i;// 


i-Pq-P^ + PqP^ - PqxP^ + HPq - PqP^ + p^ - p^p^) 


[-1 -p^-p^- p^p^ + p^p^ + i{p^ + p^p^ + p^ + p^pj^ > 


^ 27 ]L-i 2 T/L 


1 - {pQ + xp^) + i{pQ + ^c~ 


+ e^^^ [^-(1 + p^ + Pq) + i(xpQ + ^ 

1 ~{pQ + P^) + [^-(1 + </^c + ^^G + '^c + 2;//(^i//j,)sin 2 tjlJ 

+ i [(iA'G + " ^^G^c^ ■'■ ^^G 2;//gi//p)cos 27jL e^^^ + {1 + p^ + pQ)sin 2rjlj e^^^ 


III! I ■■■■ ■ II I I 


INI 



For simplicity this expression is rewritten as a + ib. Then 


-ATefa 

a + ib 

-ATqI^I + ij^^) - i;// J(cos 2r]li - i sin 2?]L)e^’^^ 

a + ib 

-AT^ + ,)cos 2 tjL - i//^ sin 2rfLl^ - cos 2 tjL + (1 + ;//^)sin 2TjLjj(a - ib) 


Rewriting the preceeding equation in simpler form using constants c and d as defined 
by equations (B30) and (B31), respectively, gives 


-ATq e^^^^(c - id)(a - ib) e^’^^ljca - db - i(da + cb)J 


A = 


27?L[ 


2 ,2 

a + b 


.2 k2 

a + b 


(B23) 


B = 


ATg[[(1 - <P^) + J ^TG(e + if)(a - ib) [ea + fb + i(fa - eb)^ 


(a + ib) 


.2 ,2 

a + b 


2 k2 

a + b 


(B24) 


where e and f are defined in equations (B32) and (B33). The solution for 6, making 
use of equations (B23) and (B24), is as follows: 


e = 


-ATg [ca - db - i(da + cb)] [ea + fb + i(fa - eb)] e“^^^ 


ICOT 


2 V.2 

a + b 


-ATq ca - db - i(da + cb)](cos tjx + i sin t/x)+ e^^[ea + fb + i(fa - eb)](cos ?]X - i sin r]x) 


a^H-b^ 
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e = ^ - db)cos tjx + (da + cb)sin jjx] + e’^^jCea + fb)cos tjx + (fa - eb)sin r;x] 


o2 u2 

a + b 


+ i I [-(ca - db)sin 77X + (da + cb)cos + j^(fa - eb)cos tjx - (ea + fb)sin Tjxje^^j^ 

The phase lag between the wall and the driving hot gas is defined as cp and is 


(B25) 


f-(ca - db)sin ?jx + (da + cb)cos rjx] + F(fa - eb)cos rpc - (ea + fb)sln nxle’^* 

(p = tan ■= L P ^ 

[(ea + fb)cos 7}x + (fa - eb)sln ? 7 xje’^’‘' - [(ca - db)cos rjx + (da + cb)sin j/xje'^^^^"^^ 


(B26) 


The amplitude ratio is 


a^+b^ 


i— '(c^' db)sin?)x +(da+ cb)cos + [(fa- eb)cos77x -(ea + fb)sin77xje’^^| 


+ |[(ea + fb)cos tjx + (fa - eb)sin T/xje*^^ -[(ca - db)cos tjx + (da + cb)sin ^ 


1/2 


(B27) 


where 


a = 1 - (i//q + i//^) + € 

+ 2i^^i//^)sin 277 L - {l + xjy^ + il^^)cos 2t/LJ 

(B28) 

i4^Q + ^c~ 

+ + 2 i//^i//^)cos 2 ? 7 L + {1 + ij/^ + ;//j^)sin 277 L] 

(B29) 


c = (1 + )cos 2rj'L - ij/ sin 2 t}L 

(B30) 


d = i// cos 2?jL + {1 + ij/ )sin 2rfL 
^ c* 

(B31) 


e=l-;^c 

(B32) 


0 

II 

(B33) 


The important equations written in nondimensional form are as follows: 
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General solution for wall temperature 


‘ Yi 




(ea+ fb)cosTjL — + (fa - eb)sin?jL — 


^ 7 ]L(x/L) 


(ca -db)cos7jL — + (da + cb)sinTjL — 


uL( 2 -x/L) 


+ 1 -! 


~ X X 

-(ca -db)sin?jL — + (da+ cb)cos?7L — 
L L 


^tjL[2-(x/L)] 


X X 

(fa - eb)cos 77L (ea + fb)sinTjL — 

L L 


5^L(x/L)|^^-io;T (234) 


or 


^ ^-i(cuT-<p) 


ATq ATq 


(B 35 ) 


where 


q)= tan 


■1 _ 


(ca - db)sin? 7 L A + (da + cb)cos? 
L 


(ea + fb)cosrjL A + (fa - eb)sinrjL A 


7 jL A - eb)cos 7 jL A - (ea + fb)sinr)L A /Taqe\ 

lJ L L y (bJODj 


rjL(x/L) _ _ db)cos 7 ;L A + (da + cb)sin 7 )L A 


j/L[ 2 -(x/L)] 


m 


‘^I'g + 




(ca-db)sin?7L (da+cb)cos tjL — 

L L 


^r?L[ 2 -(x/L)] 


(fa - eb)cosTjL — - (ea + fb)sin77L — 
L L 


,77L(x/L)|2 


(ea + fb)cos tjL — + (fa - eb)sin tjL — 
L L 


nL(x/L) 


2' 

(ca - db)cos ryL — + (da + cb)sin 77L — 

L lJ f 


, 1/2 


(B 37 ) 
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I 


Derivation of Convective Heat -Transfer Parameters as Function of 
Forcing Frequency, Wall Material Properties, and Phase Lag 

The equation which gives the convective heat -transfer parameters, and i}/^, in 
terms of the forcing frequency of the hot -gas temperature, wall material properties 
(including thermal conductivity), and phase lag angle is derived in this section. 

Equation (B26) is rewritten as 

tan (p (ea + fb)cos rjx e^^ + tan ^(fa - eb)sin rjx e^^ - tan cp(ca. - db)cos rjx 

- tan ^(da + cb)sin ?jx + (ea + fb)sin rjx e^^ - (fa - eb)cos rjx 

+ (ca - db)sin rjx - (da + cb)cos rjx = 0 

(ea + fb)(tan <p cos rjx + sin rjx)e^^ + (fa - eb)(tan cp sin rjx - cos rjx)e’^^ 

- (ca - db)(tan <p cos rjx - sin rjx)e’^^^^”^^ - (da + cb)(tan q) sin rjx + cos = 0 

(ea + fb)(sin cp cos rjx + cos <p sin rjx)e^^ + (fa - eb)(sin cp sin rjx - cos cp cos rjx)e^^ 

- (ca - db)(sin <p cos rjx - cos cp sin 

- (da + cb)(sin <p sin rjx + cos cp cos rjx)e’^^^^~^^ = 0 

(ea + fb)sin(^ + rjx) - (fa - eb)cos((jr + rjx) - [|(ca - db)sin((/r - rjx) 

+ (da + cb)cos(^ - rjx)]] = 0 (B38) 

Using equation (B38), factoring out the terms having q)^ and q>^ in them, and de- 
fining the coefficients as 

A s -sin{q> + rjx) + e^^^ cos(<p + rjx - 2rjL) + sin(^ - rjx) 

-e^^fL-x) - rjx + 2rjL) (B39) 
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B = cos{cp + rfx) - sin(^ + rjx - 277 L) - 


COS(^ - T]X) 


+ sin(^ - Tjx + 2rjlj) (B40) 

C = sin(cp + rjx) + sin((p + tjx - 2rjh) + sin((;ci - rjx) 

+ sin{(p - tjX + 2 t 7 L) (B41) 

D = -cos((p + tjx) + cos{cp + rjx - 2tjL) - cos{cp - tjx) 

+ cos(^ - TJX + 2tjL) (B42) 

E = sin((j? + tjx) - sin((;o + tjx - 2tjL) + sin(^ - tjx) 

- sin((/T - TJX + 2tjL) (B43) 

F s cos{cp + tjx) + sin (</5 + tjx - 2tjL) - cos(^ - tjx) 

- sin((jT - TJX + 2tjL) (B44) 

G s -sin(cp + tjx) - cos{cp + rjx - 2tjL) + sin(<j? - tjx) 

+ cos(^ - TJX + 2tjL) (B45) 

the desired equation is 

2(F + G)\jy^xl^^ + 2Cxjy^ + 2(C + + 2Ai//^ + (A + B)xl/^ + E - 0 (B46) 

K xj/^ is a known quantity, the equation reduces to a quadratic which is easily solved. 

K 4 /^ is known, the equation is linear and again can be solved easily. However, if 
neither 4 /q or are known, a relation must be established between them to solve the 
equation. This relation exists in the steady-state solution for heat transfer across a 
plate. 

“gK - A«»] “ - tJ <B4’) 
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Rearranging equation (B47) gives 


he _ Tq - T(0) 

ho ' T(L) - 

Defining the ratio in equation (B48) as R results in 

p.hc_ TG-T(0) 

he T(L)-T^ 

since i// = K? 7 /h 

R - — 

Substituting equation (B50) into equation (B46) gives 

o rC + R(C + D)l 9 2A + R(A + B) „ 

K + — — = — -K + - ~ ~ + — - — = 0 

^ R(F + G) ^ 2R(F + G) ^ 2R(F + G) 


(B48) 


(B49) 


(B50) 


(B51) 


The solution of this cubic equation can be obtained using the classical approach. The 
positive, nonzero, root will yield the desired convective heat -transfer parameter 
From equation (B50), then 


^G = 


(B52) 
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APPENDIX C 


DERIVATION OF THE CASE 2 SOLUTIONS - LUMPED WALL PROPERTIES 

The basic problem is the same as has been stated in appendix B. However, appen- 
dix B gives the solution for the wall temperature obtained using the distributed wall 
property differential equation. In this derivation it is assumed either that the wall pro- 
perties (density and specific heat) can be lumped and that the temperature gradient across 
the wall is negligible or that the temperature gradient across the wall (effect of finite 
thermal conductivity) can be accounted for by using a modified overall heat -transfer coef- 
ficient and lumping the wall properties at the point where the wall temperature is meas- 
ured. 

The solution for the temperature response of the wall to a sinusoidally varying fluid 
temperature is derived in the first section of this appendix assuming no wall temperature 
gradient to exist. This solution includes the exponential decaying terms (starting 
transient), an offset of the average temperatures due to a change of the average fluid tem- 
peratures, and the steady -state sinusoidal oscillations. The starting transient is of in- 
terest in determining the time required to reach steady-state oscillating conditions. 

The steady-state oscillation portion of the solution has been rearranged to give the 
heat -transfer coefficients as functions of, amor^ other things, either the steady-state 
phase lag angle or the amplitude ratio. These equations are presented in the text. 

The second section of this appendix presents the solution for the steady-state oscil- 
lating wall temperature arrived at in the first section but modified to account for the 
thermal conductivity. This is accomplished by defining an overall heat -transfer coeffi- 
cient to exist between the appropriate fluid and the point in the wall where the wall tem- 
perature is measured. The density and the specific heat are assumed to be concentrated 
at the point in the wall where the temperature is measured. The resulting solution for 
the wall temperature is then rearranged to give the convective heat-transfer parameter 
in terms of, among other parameters, either the phase lag angle the amplitude ratio. 


Response of a Convectively Heated and Cooled Wall to a Sinusoidally 
Varying Hot-Gas Temperature Neglecting the 
Effect of Thermal Conductivity 

The hot-gas temperature varies sinusoidally and can be written as 

Tq = Tq + ATq sin ooT (C 1) 
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The heat balance across the wail is formulated as 




hc(T 


Tc) 


(C2) 


Substituting equation (Cl) for Tq in equation (C2) and dividing through by pcL gives 


^ (T„ + AT^ sin CUT - T) (T - T ) (C3) 

dr pcL pcL ^ 


The solution of this first-order differential equation has been accomplished using Laplace 
transform techniques. Hence the solution yields information about the time required for 
the wall to reach a steady -state oscillating condition after the start of the hot -gas tem- 
perature oscillation. 

Making the Laplace transformation in equation (C3) yields 


h„T h„AT„ 
S^[T] -Tq=-^+-^— ^ 


cu 


pcLS 


S^ + co^ P^^ 


h„ h h T 

^[T] - — ^.^[T] + ^ ^ 


pcL 


pcLS 


(C4) 


Solving for if[T] in equation (C4) gives 


^ ^c^c , ^G^T^cu 

+ + — 


if[T] = 


pcLS 


pcL(S^ + cu^) 


S + 


h„ + h„ 
G c 

pcL 


(C5) 


Rearranging equation (C5) by dividing both the numerator and denominator by pcL gives 




V - 

— Tg + — 
pcL pcL 


‘G 


pcL 


ATqCU 


S + 


'G 


h 


(C6) 


pcL pcL, 


S S + 


G 


pcL pcL 


S + 


‘G 


pcL pcL) 


(S ) 


The last two terms in equation (C6) have been rewritten using partial fractions so that 
the inverse transforms can be found. The resulting rewritten form of equation (C6) is 
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V + h ' 

G c 




1 _ 1 \ ^ ^G^'^G 

S g^ ^G + M 
\ pcL / 


^G'‘‘*^c\ 2 A ^G‘''^c\ ^ -io^2i(S+iaj) ( ^ — + io;) 2i(S - io;) 

"" r^/ V / V / 


The inverse Laplace transform of equation (C7) yields the wall temperature 


" i'G + ''c L 


-[(h(^+hc)/pcL)]T 


-[(ho+h )/(pcL)]-) 




\ /hp + h„ \ 

^ ^ - icul 2i 1 ^ + iwi 2i 


The last two terms in the bracketed quantity in equation (C8) have been combined by re- 
writing the numerators in complex polar form and converting the resulting polar form to 
the equivalent trigonometric form as follows. 

Define an angle <p as shown in the following sketch: 




\2 -f 



hr- + h„ 
G c 
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The last two terms of equation (C8) are then written as 



Substituting this result into equation (C8) and collecting terms gives the equation for the 
wall temperature response to a sinusoidally driven hot -gas temperature with the initial 
wall temperature being T^. The response, assuming the wall material has a large 
thermal conductivity, is 


Offset from con- 

Starting transient ditions at r = 0 



tJt - ^ 




-[(h^+h )/(pcL)]r h^T„ + h T 

^ ^ ^ ~u{t) 

h^ + h 
G c 


Steady -state oscillation 
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where 


cp. stan~^ (CIO) 

>'G + *'c 

u(t) = 0 for T < 0 

u(t) = 1 for T > 0 


and the amplitude ratio at steady-state oscillation is given by 


^ml 

AT^ 


“G 

pcL 


pcL ) 


(Cll) 


+ Cl) 


Thermal Conductivity Accounted for in Lumped Solution 


It is possible to account for the thermal conductivity in the lumped solution (eq. (C9)) 
by assumii^ that the convective heat -transfer coefficient h can be replaced by a modi- 
fied overall coefficient: 




K + x^h^ 


U„ - 




^ K + x„h„ 
c c 


(C12) 


(C13) 


where is the distance into the wall from the hot -gas side and x is the distance in- 
Lr ^ 

to the wall from the coolant side. Then 


x„ = L - X, 


G 


(C14) 


The last term in equation (C9), the steady -state oscillation portion of the solution, has 
been rewritten as follows: 
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«1 = 


‘'g+'>c 


sin(coT - <pj) 


1 + 


pcLco 

^*>G + “cy 


Using equation (CIO) this expression becomes 


ATp h / 2 

^ 1 + — f 1 + tan^ cp^ 


sin(ct)T - (p-j^) 


'G 


1 

= — cos <j?H sin(ojT - <p^) 

lA 


(Cl 5) 


Substituting the modified overall coefficient U (eqs. (C12) and (C13)) for hQ and h^, 
respectively, in equation (Cl 5) yields 


cos (p2 




sin(cuT - 


'G 


K + (L ~ x^^)h , 




cos 


cp2 


ATr 


K ^G 
+ 


sin(o;T - ^2) 


1 + 


Lh^ L 


Lh, 


where subscript 2 indicates that the modified overall coefficient has been used. 
This equation has been rewritten to include ij/ as follows: Since 


jL= a 

?7L hL 


(C16) 
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and 




R 


(C17) 


09 COS cp.y 

r — sin(wT - <p^) (C18) 

x^\ 

1 

V 

This equation gives the wall temperature amplitude ratio derived from the first-order 
differential equation and accounts for the thermal conductivity of the wall material by 
using a lumped system. A lumped system in this case means that the wall is divided into 
parts and that the conductance of each part is added to that coefficient which applies to 
the boundary shared by the appropriate coefficient and the assigned portion of the wall. 

The phase lag (eq. (CIO)) has also been reworked to account for the thermal conduc- 
tivity using the lumped system approach. 

Substituting equations (C12) and (C13) for h^ and h^, respectively, in equation (CIO) 
and using the subscript 2 to indicate the resulting phase lag angle yields 

^9 = tan"^ - ^pcLco 

^ h„K h K 

c 

K + K + (L - Xj^)h^ 

Multiplying both the numerator and the denominator by 2 and rearranging this equa- 
tion gives 


AT 


G 



^2 


tan 


2h. 


pcLcu pcLco 



(C20) 
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Using equations (C16) and (Cl 7) and since 


rlyr,h=^ 

2h 


equation (C20) has been rewritten as follows; 

-1 2 


<P2 


- tan 


R 


L 





1 

0 

1 


<^2 


= tan' 


.1 V ^ ’^Gy 


1 + 


Rfl + !GnL\ 

V L 

1 7 i _! gW 

V L/ 


(C21) 


(C22) 


(C23) 


Equation (C23) gives the phase lag angle derived from the first-order differential equa- 
tion and accounts for the thermal conductivity of the wall material by using a lumped 
system. Reworking equation (C23) yields 


tan 


tan 


<P2 


q>2 


1 + 1 


V \ ^ *^0; 


|j = 2„LV/Q4 + ^°2t-') 

iJi- 

[J °\ 

V W ^G_ 


X, 




^)tjLR + + — ?jLR 

L/ L 


^ 2r,L(^Q + ^^Ll 


>^0 + 


i'-’-i) 


tjLR 




2 + 2 fl - - tan ^ 3(1 + ^) 


2x / X \ 

i // g + — ^ [1 - — jTjL^R-TjLRtancpg^O 








+ ^ (1 _ !^j,jL^R - tan (^2 = 0 (C24) 
L / 2 
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This equation relates the convective heat -transfer parameter \j/^ to 77 L, R, x/L 
and cp. Since = Ri// equation (C24) can be rewritten in terms of xj/ as 

VJ L/ 0 


+*<rjL 


1 -^ 


R 


tan (/?„ 


2 t}L 


+ /i _! gW 


R L V L/ R 2R 


From equation (C24) 


= -nk 

^G 

+ 

(1 - — )r 

2 

L 

V l/ _ 


tan <pJl + R) 1 

+ 2 ± i I 

2 • 2t/L 2 




tan ^ 2 ( 1 + 
2 qL 


0 (C25) 


1/2 


5i(l-!GWR.!itan ^2 

L \ L / 2 


(C26) 


if x^/L = 0 


OLR , ^ . 1 . 

^ 2 4tjL 


± —< 


??LR 


tan ^ 2(1 + R) 


n 2 


2r?L 




+ 2R tan 


^2} 


1/2 


if Xq^/L = 1 


»^G = 


_^L ^ 2(1 + R) 1 


4 rjL 


± -■< 

2 


i]L - 


tan ^ 2(1 + R) 


v~i2 


2 rjL 


+ 2R tan ^2 


1/2 


The amplitude ratio taken from equation (C18) has been rearranged to give as 
follows: 


^m 2 


AT 


G 


r 

/ ,, \ 

~i 

/ V 


r' 



1 

1 . ( 1 -::°) 

tjLR 

^.R fl + !G3k\ 

*■ - 

1 + 


tjLR 

1 

\ w 

^G_ 

V L 4^J\ 
J 


- 

V L/ 

*^G_ 


cos f /?2 = 0 


*^m 2 

AT^ 


X 


G 


X 


G 




G 


X 


+ II - t)LR 


cos (p2= 0 
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This equation gives the convective heat -transfer coefficient i//^ calculated using a 
lumped wall property system. 
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